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Abstract 

Let r{{p) denote the orbit of Mercury. We compare the formulae obtained via general 
relativity for r{p) and for the corresponding perihelion precession angle Ai^, with the 
formulae obtained via the relativistic gravitational law, F = GMm'y^/r'^. The latter 
law can be derived from Newton’s gravitational law by employing the gravitational 
rather than the rest masses of the Sun and Mercury. Remarkably, it is found that the 
two expressions for the dominant part of r{ip) are identical and that the two expressions 
for A(f are also identical. 
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1 Introduction 

The experimental verification of the formulae for the perihelion precession of Mercury and 
for the bending angle of light passing near massive bodies such as our Sun, are the two 
most spectacular confirmations of Einstein’s theory of general relativity [T]. Furthermore, 
the efforts of Einstein in order to explain these two phenomena played an important role 
in the development of general relativity (GR) [2]. It is generally accepted that Newton’s 
universal gravitational law, which cannot account for the observed perihelion advancement 
and which gives only half of the experimentally observed light bending angle, is not accurate 
enough under relativistic conditions nn. 

However, when employing Newton’s universal gravitational law under relativistic con¬ 
ditions, little attention has been paid to the fact that there exist differences between the 
rest, relativistic, inertial and gravitational masses, of the bodies involved. 

The inertial mass, m^, and the gravitational mass, nig, are equal according to the equiva¬ 
lence principle [5] , but the ratio of the inertial mass to the rest mass and also the ratio of the 
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relativistic mass, 'yruo, to the rest mass are unbounded as the particle velocity approaches 
the speed of light [6] . 

Consequently, the following important question arises when using Newton’s universal 
gravitational law under conditions where the Lorentz factor 7 is significantly larger than 
unity: should one use the particle rest mass, nio, the relativistic mass, yrUo, or the gravi¬ 
tational mass, nig, which equals the inertial mass, ruj? This question has been addressed 
recently in [3 El [g do], where it is argued that one should use the gravitational mass, nig, 
rather than the rest mass, ruo, or the relativistic mass, yrUo. Perhaps one of the reasons that 
this question was not discussed earlier is that the velocities of the planets of our planetary 
system are of the order 10 “^c, thus the corresponding Lorentz factors are of the order of 
1 -|- 10“® and therefore the differences between rest, relativistic, and gravitational masses 
are very small, which is also the case for the magnitude of the perihelion precession value 
per orbit and of the light bending angle caused by our Sun. However,this is not the case 
for many other astronomical systems involving heavier stars mm- 

The inertial mass, m*, and hence according to the equivalence principle, the gravitational 
mass nig, is related with the rest mass, mo, via the equation 

mg = mi = -f^mo, 7 = (1 - (1) 


This result for linear motion was first obtained in Einstein’s pioneering special relativity 
paper IS nail]; it has been argued recently mu that this result remains valid for an 
arbitrary particle motion, including circular motion. 

Consequently, Newton’s universal gravitational law for the gravitational force between 
two particles, at a distance r, of rest masses mo,2 and moving with a velocity v relative 
to each other, is given by 

^ _ Gmg^img^2 _ Gmo,imo,2'y^ 

^ 2 2 ' 


If mo7 = mofi = mo, the above equation reduces to 

Gml-f^ 


( 3 ) 


By employing equation Q it is possible to predict the existence of gravitationally con¬ 
fined circular rotational states with many of the properties of hadrons mum, or bosons 
m, and to show that the rest masses of quarks are very small, in the mass range of neutri¬ 
nos mM- Furthremore, using equation ([ 3 ]) it can be shown that the ratio F/Fo of the force, 
F, keeping two or three-particle systems in orbit, divided by the same force computed for 
7 = 1 , is {mpi/mo)^, which is exactly the value predicted by GR [U [8]. 

For the motion of a particle with rest mass mo rotating around a massive particle of 
rest mass Mo, equation (jj]) remains valid. Indeed, consider for example the Sun (mass Mo) 
and Mercury (mass mo) at the locations of two instantaneous inertial systems, S and S', 
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Figure 1: Model geometry showing the coordinate system and the definitions of (po and p. 

(Fig. 1) moving with an instantaneous velocity v relative to each other. In this case, the 
inertial and thus the gravitational inertial mass of Mercury is given by 'y^rrio according to 
the observer on the Sun at S, but also the inertial mass of Sun is given by according 

to the observer on Mercury at S'. 

If one uses rest masses instead of inertial masses then, according to Newtonian mechanics 
the particle of mass rrio moves on an ellipse. Consider such an ellipse with an elliptical 
parameter p (also known as the semi-latus rectum), with eccentricity e, and with a semi¬ 
major axis a; these parameters are related by the equation 

p = {l- e‘^)a. (4) 


GR implies the following expression for the precession angle, Ap, [I5lll6l[171[l8l[l9l[20] 


Ap 


evrGM 

pc^ ’ 


J2 

GM’ 


(5) 


where J denotes angular momentum per unit mass. 

Using the momentum-velocity equation of special relativity (SR), namely momentum 
equals ynioV, Lemmon and Mondragon examined in [21] the resulting corrections to Keple- 
rian orbits and to the precession angle Aip. Their work, which does not take into consider¬ 
ation equation m, yields one sixth of the value computed via GR (equation Q). 

It is interesting to note that Silberstein already back in 1917 [22| attempted to account 
for the velocity dependence of the rotating mass m by combining m = ymo with the 









empirical force expression 


GMuio ,, . 
F= -o— fiV- 


( 6 ) 


He reported that he could get agreement with the GR result of equation ([5]) provided 
that 7 ( 7 ) = 7 *^“^ and n = 6. He noted that it was entirely unclear to him why this 
particular value of n leads to exact agreement with GR |22j . Furthermore, Silberstein only 
analysed the relevant equations in the limit of v/c —)• 0, where v is the velocity of the 
rotating particle. 

Here, for the sake of completness we first present a simple rederivation of the well known 
d la isi ttzi dsi da [20] GR expression 


d‘U , /3e 




(7) 


where 

rr j L 2GM r, ^ (GM\ 

U = 77 ^ , J = — , rs = —— , e = — = 2[ —- 
GMr m & p \ cJ J 


( 8 ) 


In equations dSj)! R is the angular momentum of the particle of mass m, M is the mass of 
the heavier body, G is the gravitational constant, and J is the angular momentum per unit 
mass. 

We then present the first main result of our paper: if a particle of mass m moves around 
a heavier mass M under the influence of the relativistic gravitational law ([2|) , then U satisfies 
the ordinary differential equation 


(fU 

d(p‘^ 


+ U 




(9) 


where _ _ _ 

r~_ <72 L ~_rs_^(GMy 

^ GM? ’ mo' ^ p V cJ J ’ ^ GM 

and ~ denotes quantities computed via equations l|T |) to ([3 ]) . 

It turns out that for small s the dominant parts of the solutions of equations ([7|) and 
(d are both given by the expression 


Ud 



etp sin 7 , 


( 10 ) 


and thus, the precession angle per revolution, Aip, is given in both cases by 

= 2^ (I) . 


( 11 ) 


which using the definition of e in the third of equations (d) becomes the well known formulae 













It is important to note that one may consider a third instantaneous inertial system, S", 
with an observer, located on earth and with a velocity v/j with respect to one of the two 
moving bodies. This second velocity v^j, and the corresponding Lorentz factor, yjj, dehnes 
the ratio 

nir/mo = 'jR, ( 12 ) 

of the relativistic and rest masses of Mercury as observed from S", i.e. from the Earth. 
Since nio is not amenable to direct observation, and also since Aip should be independent 
of S", we may set jr = 1. 


2 General Relativistic Treatment Revisited 


Let (r, (p) be cylindrical coordinates in a plane containing the center of the sun, and let r, 
E, and L denote the proper time, energy, and angular momentum respectively. 

The basic equations of GR are conservation of energy and of angular momentum, as 
well as the Schwarzschild metric equation: 


E 

mc^ 



dt 

dr' 


(13) 


L 


9 dip 

m dr 


(14) 


c^dr^ = cM 1 - 


dr- l- 


dr^ — r'^dp^, 


2GM 


r. = 


(15) 


(16) 


where is the Schwarzschild radius of the sun and M is the mass of the sun. 

Solving equations (fTSl) and (fTTl) for dt and dr in terms of dp, substituting the resulting 
expressions into (1151) . and then solving the resulting equation for (dr/dp)"^ we find 


I dr E"^ / 

dp J 


~1X 


1 

H —^ 


(17) 


Making the substitution 

equation (fT7)l becomes 


du 

dp 


u = 


E^ 


- (1 - TsU) 


~1X 


+ r 


(18) 


Differentiating this equation with respect to p and introducing the dimensionless quan¬ 
tities dehned in equations ([8]) we find equation ([7|). 








3 Motion Due to a Central Force 


Consider a single particle of relativistic mass and of gravitational mass rUg moving under 
the influence of a central force of magnitude /. Then, 


— (m^f) = -/r, 



(19) 


Incidentaly, in the particular case that nir = mo"), using the definition of 7 , i.e. the 
second of equations m, we find 


dj 7 ^ dv 
dt dt ’ 


V = If 


( 20 ) 


Hence, 


dt 


{mo'jr) 


mo'yr + 


mo")^ 



( 21 ) 


It is shown in the Appendix that 


dv 

V— = r ■ r, 

dt 


and then equation () 2 ip takes the familiar form found in most books, namely 


d TTlo")^ 

— (mo'yr) = mo'yr H- 5 — (f • f) f • 

dt 


( 22 ) 


Starting with (I19p . it is straightforward to derive the following results: 


1. The angular momentum L is conserved, where L is dehned by 


L = rx mrr. 


(23) 


Furthermore, the magnitude of the angular momentun denoted by L is given by 

L = (24) 

where (r, (p) are cylindrical coordinates in the plane of motion which clearly is orthog¬ 
onal to L. 

2 . u{ip) satisfies the equation 


d^tt fmr 



3. The speed of the particle v satisfies the equation 


v 


2 



(25) 


(26) 







In order to derive the above results we begin by differentiating equation (|23p with respect 


to t: 


L = rx nirZ + r x -j^{mrr) = r x {—fr) = 0. 


(27) 


Thus, L is constant. 

Equation (I23p implies that r is perpendicular to L. 

Introducting cylindrical coordinates in the plane perpendicular to L, we find 


r = rr, f = cos (pi + sin pj, (28) 

where i and j are unit vectors along the x and y axis. 

Differentiating (1281) with respect to t we find 


r = {r cos p — rep sin p)i + (r sin p + rp cos p)j. 


(29) 


Thus, the definition of L implies 


L = m. 


•r‘^(pK, 


(30) 


where k is the unit vector perpendicular to L. 

Let L denote the magnitude of L. Equation (1301) implies 


L 

p = -; 

nirV 


(31) 


Decomposing equation ([TUp along the x and y axis we find the following equations; 


^{mrx) = -fcosp, ^{niry) =-f sin p. 


(32) 


Using the indentity r = 1/u we find 


nirX = mr — {r cos p) = mrp—{u ^cos;^). 
at dp 


(33) 


Replacing in this equation p via equation (|3ip and simplifying we find 


rrirX = —L 


( du\ 

usinp + cosoj— 

V dpj 


( 34 ) 


Thus, the first of equations (l32l) yields 


—Lp 


du . . du d?u 

u cos p + —s\np — sm p- -h cos p-r^ 

dp dp dp^ 


—f cos p. 


Replacing in this equation (p by m, and simplifying we find equation (I25p . 




(35) 


Equation (I29p implies that v = |r| satisfies 


V 


2 


•2 , 2-2 
r +r ip . 


Noting that 


dr d _i _2 du 

r = — = (p—u =-(pu 

dt dip dip 


and using for ipu ^ the expression obtained from (|3ip . namely 


ipu 


L 

nir ’ 


it follows that 

L du 
rur dip 

Using in (l35P equations p6P and (f37P we find eq. ([261) . 


(36) 


(37) 


4 The perihelion precession of the Mercury 

In this case 

rrir = irio, f = Grtio'^^u^ ■ 

Thus, equations (f25P and ([261) become 

d’^u _ GM fj 

+ u = -^7^ 


dp^ J2 


and 


v 2 = P 


■••'ST 


where J denotes angular momentum per unit mass, i.e. 

L 


J = 


rUn 


Introducing the dimensionless parameter U via 


•7 ^ p 


equations (1391) and (|10P become 


fu_ 

dp'^ 


+ t/ = 7® 


and 


V _ e 

7 ^ ~ 2 


U^ + 


dp 


(38) 


(39) 


(40) 


(41) 


(42) 


(43) 


(44) 










where £ is defined by 


?=^ = 2 f™V. 

P \ cJ J 


(45) 


Equations (1431) and (1441) together with the definition of 7 imply 

-3 


~ I , e 

I 2 


In order to solve (1461) for small e we let 




dip 


(46) 


U — Uq + ^(^)) ^ ^ 0- 


Then, equation (l46]) yields 


'^ + U =1 
< i »)2 ° ’ 


and 




dp'^ 

The general solution of equation (|48p is 

Uo = 1 + ecos{p - po), 

where e and po are constants. Letting for simplicity po = 0, equation (14^ becomes 

d?Ui 


(47) 

(48) 

(49) 

(50) 


dp'^ 


+ C/i = 1 + e + 2e cos p. 


Thus, 


Ui = 1 + + ep sin p. 

Hence, using (l50l) with (/?o = 0 and (15^ into (1471) we find 

~ 3 

U = 1 + e cos p + -£\l + + ep sin p\ . 

The dominant part of U denoted hy Ud is given by 


Ud = -ee(/9sin p. 

We next solve equation ([7|) for small e. Letting 


(51) 

(52) 

(53) 


(54) 


U = Uo + -£Ui, 


( 55 ) 










we find that Ui satisfies the ordinary differential equation 


d^Ui 2 2 

, „ + 17i = 1 + e cos ip + 2e cos ip. 
dip^ 


Thus, 


U = 1 + e cos ip + -e 


3 3 

The dominant part of U denoted by Ud is given by 


1 H—-cos^ ip + eip sin ip 


( 56 ) 


(57) 


Ud = -eipsm.ip. 


(58) 


In order to relate e and e we need to relate J and J. In this respect we note that 
equations (fTTl) and (1^ yield 


T = r^^ 
dr ’ 


7 ~>2dip 
J = r 

dt 


(59) 


Thus, 

J (^ dt 
J \u) dr' 

Hence, using equation (fT^t) to compute dt/dr, as well as using the first of equations ([S]) and 
equation ()42p to relate u and u with U and U we find the following: 


J _ (uV f jV E/mc^ 

j~[u) T^’ 


or 


E U 


mc^ \ U j 1 — eU 


(60) 


Equation (IlSp can be rewritten in the form 


dUy _2 f E 
dip) 


e \ mc^ 


-il-£U){- + U' 


(61) 


Using in this equation the expansion 


U = 1 + ecos ip + 0(e), 


it follows that 

—^ = 1 + 7(e^ - 1) + O(e^), e 0. 
mc2 4 

Letting in (l60P 

E 

—2 = 1 + 0(e), £ —>■ 0, 

mc^ 


( 62 ) 












Figure 2: Plots of the numerical integration of equations ([7]) (general relativity, GR) and 
([9]) (relativistic gravitational law, RG) and comparison with the Keplerian orbit (j48jl for 
reference. Distances in 10^® m. For illustrative purposes the value used for GM is larger by 
a factor of 10^ than the real value of 1.3271 x 10^*^ rr? j. 


and 


we find 


Hence, 


— — 1 + 0(6',?), E —>■ 0, E —> 0, 

= 1 - 1 - 0 ( 6 , 6 —y 0 , 6 —^ 0 . 

J 

E = e(l -l- 0(6)), 


(63) 

(64) 


thus Ud coincides with Ud- 

The computation of the precession of the perihelion depends only on Ud- Indeed, 
equation ([57)1 implies that dU/dtp vanishes at (^ = (^*, where satisfies the equation 


sin P* = [—-E sin 2^^, -|- sin (/?* + cos p* ) . 


(65) 


One solution of this equation is = 0. In order to obtain the second solution we let 
Pif = TT + Eu and we use the identities 


cos(7r-|-60;) = —1-|-0(6^) , sin(7r-|-6o;) 


EUJ + 0(6^), 6^0. 


( 66 ) 












Then, equation (l65]) implies 


3 3 TTTs 

¥?* = vr + -7re = tt + -- = vr + 

2 2 p 


2,'kGM 

pc^ 


(67) 


Equation (f53]i with e = s also yields for (p^ an expression identical with (f 66 l) . 

Results of numerical integration of equations ([7]) and Q as well as of the Kepler equation 
(f78l) are presented in Figure 2, which shows that indeed the results of equations ([7]) and 
([U|) coincide for small e, i.e. the 7 ® gravitational law yields to order e the same results with 
general relativity. 


5 Conclusions 

The present results show that the relativistic Newton’s gravitational Law, which uses grav¬ 
itational rather than rest masses and thus contains the term 7 ^ multiplying each moving 
mass in the familiar Gmim 2 /r‘^ or GMm/r'^ term, leads to a particle trajectory r{ip) whose 
dominant term is the same with the dominant term of the familiar r{ip) expression obtained 
via the theory of General Relativity. Furthermore, the perihelion precession angle, Aip, 
computed via the former approach is equal to the value computed via General Relativity. 
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APPENDIX A 


Here we establish the validity of the equation 

dv . .. 

V— = r ■ r. (A-1) 

dt - ^ ^ 

Differentiating equation ([2^ with respect to t we hnd 

r = rr + (2r0 + r(p)r^ — r(p‘^r = (r — r0^)r + (2r<^ + r(p)f^, (^’2) 

where the unit vector f"*' define by 

— sin ifi + cos (pj, (A-3) 

is orthogonal to the unit vector r. 

Equations (l28l) , (f29]l , and (|A-3I) imply that r can be rewritten in the form 

r = rr + rpr^. (A-4) 

Hence using r-f^ = 0, it follows that 

r • f = (r — r(p‘^)r + {2rip + r(p)rip. (A-5) 


Thus, 


. .. ... ■ ■ 2 2 • 

r • r = rr + rrp + r pip. 


Differentiating (1351) with respect to t we hnd 


dv 2 2 

V— = rr + rfp + r pip. 
dt 


(A-6) 


(A-7) 


Equations (IA-61) and (IA-7P imply (lA-ip . 










APPENDIX B 


We will establish the following fact: suppose that r satishes the equation 


Then, H dehned by 


= - ft 


n , , dv 

= "io 7 C + V(r), /=^> 


is constant. 

Indeed, equation (|B-2p implies 


dH 

dt 


, dq dV dr \ 
dip dr dip J 


—— = ( moc^-t- + -r~r ) ‘f- 


Recall that 


Thus, 


2 . 

7 =1 + -T 


d7 




dip J 


u + 


d 2 


u 


du 

dip 


Furthermore, u satishes 


d2 


U 


dip" 


+ u = 


dip"^ 

7 dV 
mo'Pu^ du 


Hence, equation (IB-41) becomes 

dH rrioC^ 7 dV dV d _i 
dt [ 7 moHu"^ dr dr dip 


ip = 0. 


(B-1) 

(B-2) 


(B-3) 


(B-4) 
















